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Scaling behavior and positivity violation of the gluon propagator in full QCD 
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The Landau-gauge gluon propagator is studied using the coarse and fine dynamical MILC con- 
figurations. The effects of dynamical quarks are clearly visible and lead to a reduction of the 
nonperturbative infrared enhancement relative to the quenched case. Lattice spacing effects are 
studied and found to be small. The gluon spectral function is shown to clearly violate positivity in 
both quenched and full QCD. 
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I. INTRODUCTION 

Quarks and gluons, which carry color charge, are not 
observed as free particles, but only inside the colorless 
bound states called hadrons. This phenomenon is called 
confinement. The gluon propagator has long been stud- 
ied regarding confinement. It is also an important in- 
put for many phenomenological calculations in hadronic 
physics (see e.g. the reviews [1-6]). 

In the last decade there have been many lattice studies 
devoted to the gluon propagator in the Landau gauge. 
Most of them are done either in quenched QCD [7-20] 
or in quenched SU(2) [21-23], but there are also more 
recent reports [24-27] by us and others for the case of 
full QCD. All lattice studies so far have indicated that 
the gluon propagator in Landau gauge is infrared fi- 
nite (see e.g. [14, 24, 28]). This, however, disagrees 
with the corresponding results obtained in studies of the 
Dyson-Schwinger equation (DSE) for the gluon propa- 
gator. There an infrared vanishing gluon propagator is 
predicted [1, 29-31]. Such an infrared behavior has also 
been found using stochastic quantization [32] and in stud- 
ies of exact renormalization group equations [33] . Recent 
DSE studies on a finite torus strongly suggest that this 
discrepancy is due to finite volume effects [34] . 

In this paper we extend the previous work reported 
in [24] to finer lattices studying the scaling behavior of 
the gluon propagator in Landau gauge. We show that 
the gluon Schwinger function in quenched and full QCD 
is negative in a certain interval. This is consistent with 
both a recent lattice study of two flavor QCD [28] and 
with results obtained in DSE studies [1, 35]. 

We use a one-loop Symanzik improved gauge ac- 
tion [36]. It was noted some time ago (e.g. Ref. [37]) that 
improved gauge actions can lead to violations of positiv- 
ity. These violations are lattice artifacts and vanish in the 
continuum limit. In contrast, the positivity violations of 



the gluon propagator reported here are independent of 
the lattice spacing and are therefore expected to survive 
the continuum limit. They are not associated with the 
choice of action. 



II. POSITIVITY OF EUCLIDEAN GREEN'S 
FUNCTIONS 

Here, a note on positivity and its violation is in order. 
In quantum field theory in Minkowski space, if a cer- 
tain degree of freedom is supposed to describe a physical 
asymptotic state, it must not have any negative norm 
contributions in its propagator. That is, the propagator 
must not violate positivity. Otherwise the states it de- 
scribes cannot be part of the physical state space; they 
are, so to say, confined from the physical world. 

Considering Euclidean Green's functions, positivity 
translates into the notion of reflection positivity as one of 
the famous Osterwalder-Schrader axioms [38, 39] of Eu- 
clidean quantum field theory (see also [40, 41]). For our 
purpose it is instructive, and also sufficient, to consider 
reflection positivity in the case of an Euclidean 2-point 
function, D(x — y), whose corresponding propagator in 
momentum space, D(q 2 ), can be written in a spectral 
representation 



D{q 2 ) = J dm 2 



p{m 2 ) 
q 2 + m 2 



(1) 



where q 2 denotes the four-momentum squared in Eu- 
clidean space and m 2 is the mass squared. The spectral 
function, p(m 2 ), is unknown in general, but if p(m 2 ) > 
for all m 2 , Eq. (1) is known as the Kallen-Lehmann rep- 
resentation. 

The 1-dimensional Fourier transform of D(q 2 ) at zero 
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spatial momentum defines the wall-to-wall correlator [42] 

oc 

C{t) = j dmp(m 2 )e- mt . (2) 
o 

In lattice spectroscopy the exponential decay of this func- 
tion is used to extract the mass of a particle state. 

Obviously, from Eq. (2), if the spectral function is pos- 
itive, then C(t) > 0. This needs not to be the case the 
other way around. If, however, C{t) is found to be nega- 
tive, p(m 2 ) cannot be a positive spectral function. That 
is, there is no Kallcn-Lehmann representation and the 
corresponding states cannot appear in the physical par- 
ticle spectrum: they are confined. 

Alternatively, an effective mass [7, 42] 

m cff (t) := --|lnC(i) 

could be considered. For positive p, the slope of m e s(t) 
cannot be positive. If, on the contrary, m e s(t) increases 
with t, p cannot be a positive spectral function. 

The gluon 2-point function was known from the begin- 
ning to violate positivity. Already in the first numeri- 
cal study [7] of the gluon propagator in quenched lattice 
QCD in Landau gauge this has been demonstrated by 
an effective gluon mass rising with increasing distance t 
(see also [43-45]). A direct observation, however, of a 
non-positive gluon correlator, C(t), in particular for full 
QCD, has been unfeasible for a long time. Nevertheless, 
in more recent times a non-positive C(t) has been ob- 
served in studies of the corresponding Dyson-Schwingcr 
equation [35], in quenched SU(2) lattice gauge theory in 
three dimensions [46] and hinted at in QCD with dynam- 
ical quarks [47]. Explicit evidence in quenched and full 
QCD (clover-improved Wilson fermions) has been given 
recently in [28]. To complete those findings, in this paper 
we show that C(t) < for some t in full QCD with 2+1 
flavors of 0(a 2 ) Symanzik-improved staggered fermions. 

III. GLUON PROPAGATOR ON THE LATTICE 

In the continuum, the Euclidean gluon propagator in 
Landau gauge has the tensor structure 

DfM={s^- q -f-y^D( q 2 ). (3) 

Here D{q 2 ) is a scalar function which contains the whole 
nonperturbative information of D^l(q). At tree-level 

D(l 2 ) = ±. (4) 

Assuming Eq. (3) to be well satisfied on the lattice, 
we have calculated D(q 2 ) using a variety of quenched 
and dynamical configurations. These configurations were 



TABLE I: Lattice parameters used in this study. The dy- 
namical configurations each have two degenerate light quarks 
(up/down) and a heavier quark (strange). In the table we 
show the sea quark masses both in dimensionless (lattice) 
units and estimated physical units. For details see Ref. [48]. 





Dimensions 


P 


a (fm) 


ma m (MeV) #Config. 


1 


28 3 


x 96 


8.40 


0.086 


— quenched — 


150 


2 


28 3 


x 96 


7.09 


0.086 


0.062, 0.031 14, 68 


108 


3 


28 3 


x 96 


7.11 


0.086 


0.124, 0.031 27, 68 


110 


4 


20 3 


x 64 


8.00 


0.120 


— quenched — 


192 


5 


20 3 


x 64 


6.76 


0.121 


0.010, 0.050 16, 82 


203 


6 


20 3 


x 64 


6.79 


0.120 


0.020, 0.050 33, 82 


249 


7 


20 3 


x 64 


6.81 


0.120 


0.030, 0.050 49, 82 


268 


8 


20 3 


x 64 


6.83 


0.119 


0.040, 0.050 66, 83 


318 



provided to us by the MILC collaboration [48] through 
the Gauge Connection [49] and were generated with the 
0(a 2 ) one- loop Symanzik improved gauge action [36]. 
For the dynamical configurations the "AsqTad" quark 
action was used. This is an 0(a 2 ) Symanzik-improved 
staggered fermion action with 2+1 flavors implemented 
using the "fourth-root trick". The lattice spacing was 
determined from the variant Sommcr parameter, n [48]. 
Following MILC convention we refer to the 28 3 x 96 con- 
figurations as being "fine" lattices and the 20 3 x 64 config- 
urations as the "coarse" lattices. See Table I for details. 

With this lattice gauge action the Landau-gauge gluon 
propagator at tree-level is 

m j 

where 

^ = ^f' n " e {-2>-2\> (6) 

a is the lattice spacing and is the length of the lat- 
tice in the p direction. To ensure the correct tree-level 
behavior for the lattice gluon propagator, a "kinematic" 
choice of momentum, 

is employed [14]. 

Lattice Monte Carlo estimates for the bare gluon prop- 
agator, D(q 2 ), have to be renormalized. The renormal- 
ized propagator Dn(q 2 ;p 2 ) is related to D(q 2 ) through 

D(q 2 ) = Z 3 (p 2 ,a)D R (q 2 ;p 2 ) (8) 

where p is the renormalization point. The renormaliza- 
tion constant Z 3 depends on the renormalization pre- 
scription. We choose the momentum space subtraction 
(MOM) scheme where Z 3 (p, 2 ,a) is determined by the 
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FIG. 1: Gluon dressing function in Landau gauge for the fine 
lattices (28 3 x 96). Full triangles correspond to the quenched 
calculation, while open circles correspond to 2+1 flavor QCD. 
The bare light quark mass is m = 14 MeV for the full QCD 
result. As the lattice spacing and volume are the same, the 
difference between the two results is entirely due to the pres- 
ence of quark loops. The renormalization point is at /i = 4 
GeV. Data has been cylinder cut [16]. 



tree-level value of the gluon propagator at the renormal- 
ization point, i.e. 



D R (q 2 ) 



1 



q 2 =fl 2 A 1 



(9) 



In this study, we will choose to renormalize either at 
(i = 1 or at [i = 4 GeV. 



IV. RESULTS 
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FIG. 2: The sea quark mass dependence of the Landau gauge 
gluon dressing function for the fine lattices. Filled triangles 
correspond to bare light quark mass m = 27 MeV, open circles 
correspond to the bare light quark mass m = 14 MeV. Data 
has been cylinder cut [16]. No mass dependence is observed 
for this case. 



the lightest). It would be interesting to study the gluon 
propagator with heavier sea quarks so that the transition 
between quenched and full QCD could be observed. We 
leave this for a future paper. 



B. The scaling behavior 

The renormalized propagator becomes independent of 
the lattice spacing as the continuum limit is approached. 
We can see this expected result in the renormalized 
quenched dressing function in Fig. 3. There data from 
the coarse and fine lattices are compared, having been 



A. The effects of dynamical sea-quarks 

To begin with, we discuss the effects of dynamical 
sea-quarks on the gluon propagator. First we compare 
the dressing function, q 2 D(q 2 ) 1 in full QCD to that 
in quenched QCD. For the fine lattices this is shown 
in Fig. 1. Obviously, the infrared hump seen in the 
quenched case is somewhat suppressed (about 30%) due 
to color screening by the quark-anti-quark pairs. The 
basic shape, however, is the same in the quenched and 
dynamical cases. This is consistent with previous re- 
sults on smaller, coarser lattices [24] by some of us, and 
also with results obtained independently using clover- 
improved Wilson fcrmions [25, 26]. 

In Fig. 2 we are looking for the dependence of the dress- 
ing function on the sea quark mass. In this case the 
light quark masses differ by a factor of two (the strange 
quark has the same mass in both cases), but we see no 
effect. Note that in Ref. [24] a small quantitative dif- 
ference has been observed by studying a slightly greater 
range of masses (the heaviest had four times the mass of 
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FIG. 3: The scaling behavior of the quenched gluon dressing 
function renormalized at fi = 4 GeV. Triangles corresponds 
to the gluon propagator on the fine ensemble, open circles to 
the coarse ensemble. Good scaling is observed. 
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TABLE II: Quenched configurations used in Ref. [14] listed in 
order of increasing volume. For this study the lattice spacing 
has been set using the string tension ^fa = 440 MeV. 





Dimensions 


13 


a (fm) 


Action 


#Config. 


1 


16 3 


x 32 


4.38 


0.165 


Tree-level Symanzik 


100 


2 


32 3 


x 64 


6.00 


0.0982 


Wilson 


75 


3 


12 3 


x 24 


4.10 


0.270 


Tree-level Symanzik 


100 


4 


10 3 


x 20 


3.92 


0.35 


Tree-level Symanzik 


100 


5 


16 3 


x 32 


3.92 


0.35 


Tree-level Symanzik 


100 



0.0 



3 4 5 
q (GeV) 



FIG. 4: The scaling behavior of the renormalized propaga- 
tor in full QCD at fi = 4 GeV. Triangles corresponds to 
gluon dressing function from the fine ensemble for sea quark 
mass m — 14 MeV. The open circles is for sea quark mass 
m — 11 MeV which is obtained by extrapolating the coarse 
ensemble data. Some difference is seen in the large momen- 
tum region. 



renormalized at /j, = 4 GeV. These two sets of data al- 
most lie on the same curve, and hence we conclude that 
good scaling is found for the quenched results. 

Turning now to full QCD, in Fig. 4 we compare the 
gluon dressing function from the fine ensemble for the 
lightest available sea quark mass (14 MeV) with that at 
an approximately equivalent quark mass from the coarse 
ensemble. Ideally while altering the lattice spacing we 
should hold the running quark mass (rather than bare 
quark mass) constant. The quark mass function was 
studied in Ref. [50] using these gauge configurations and 
it was found that a bare quark mass of 14 MeV on the 
fine ensemble had the same running mass as a bare quark 
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FIG. 5: The quenched gluon dressing function renormalized 
at 1 GeV using a variety of lattice spacings and gauge actions. 



of 11 MeV on the coarse ensemble. We therefore per- 
form the small extrapolation of the coarse ensemble gluon 
propagator to 11 MeV for this comparison. Some very 
small systematic differences at the highest momenta are 
suggested, but this is where the discretization errors are 
expected to be large. 

In the quenched case we can compare with a wider 
range of data sets. In Fig. 5 we include results 
from Ref. [14]. These were produced using the 0(a 2 ) 
Symanzik tree-level improved gauge action, except for 
the (3 = 6.0 configurations which were generated with 
the single-plaquette Wilson gauge action. We renormal- 
ize at 1 GeV so as to accommodate the coarsest lattices. 
We use the string tension, y/a = 440 MeV, to set the 
lattice spacing. For (3 = 6.0 the precise measurement of 
Ref. [51] was used. A summary of these details is pro- 
vided in Table. II. The three gauge actions approach 
the continuum limit slightly differently: for the tree-level 
improved gauge action the ultraviolet tail of the gluon 
propagator drops as the lattice spacing shrinks while for 
the one-loop improved action the tail rises. They do ap- 
pear to be converging to the same result. 

We also revisit the deep infrared region of the gluon 
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FIG. 6: The quenched gluon propagator renormalized at 
1 GeV using a variety of lattice spacings and three differ- 
ent gauge actions. The value at zero four-momentum drops 
as the lattice volume increases. 
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propagator itself in Fig. 6. The 16 3 x 32, = 3.92 lattice 
has by far the largest physical volume (4.93 3 x 9.86 fm ) 
and the smallest value for D(q 2 = 0). The propagator 
systematically increases as the four- volume gets smaller. 
Finite volume effects are no longer significant for mo- 
menta above about 700 MeV. 



C. Violation of reflection positivity 

Now we turn to the Euclidean time correlator, Eq. (2). 
On the lattice this function can be evaluated using the 
discrete Fourier transform 

L 4 -l 

C(t) = -= e- 2mit ' Ll D(p A , 0) , (10) 

V 1,4 «4=0 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 



quenched-28 3 x96 
m-27.0 MeV-SB'xge 
m-14.0 MeV-38'x96 



2 3 
t (fm) 



m = 16 MeV-20 3 x64 
m = 14 MeV-28 3 x96 



FIG. 7: The real space propagator C(i) for 28 3 x 96 lattices 
plotted as a function of dimensionful t for both light sea quark 
masses and the quenched case. Reflection positivity is clearly 
violated in all three cases. 
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FIG. 8: The real space propagator C(t) for the coarse and 
fine lattices at approximately same light sea quark masses. 



where L4 is the number of lattice points in time direc- 
tion, P4 is the Euclidean time component of the lattice 
momentum, 77,4 is an integer and D{pi, 0) is the gluon 
propagator in momentum space at zero spatial momen- 
tum. 

In Fig. 7, C(t) is shown for the fine lattices in the 
quenched case and for both choices of sea quark. In 
all cases the Schwinger function clearly becomes neg- 
ative signaling explicit reflection-positivity violation by 
the gluon propagator. 

As mentioned in the introduction, some violations of 
positivity are a regrettable consequence of the altered 
cut-off effects of improved actions [37]. These effects, 
being ultraviolet effects, can cause violations of positivity 
in two-point functions at short distances, particularly if 
simulations are performed too far from the continuum 



100 



N f = 2+1, rn = 14 
m = 27 
N f = 0, /3 = 8.4 




10 15 

t [1/GeV] 

FIG. 9: The top figure corresponds to the absolute value of 
the gluon Schwinger function from our fine lattice calculations 
in quenched QCD as well as in full QCD. Error bars are not 
shown for simplicity. The bare quark masses for full QCD 
simulations are m = 27 MeV and m = 14 MeV, respectively. 
The bottom figure is taken from [35]. It shows the numer- 
ical results for the absolute value of the Schwinger function 
from the DSE result compared to the fits in the infrared. See 
Ref. [351 for more details. 
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limit. Our study concerns the long distance behavior of 
the gluon two-point function, and many preceding studies 
(e.g. [14, 18]) have indicated that the infrared behavior 
of the Landau-gauge gluon propagator is robust against 
changes of the lattice spacing, at least for the range of 
parameters we use. Therefore, our results are unlikely 
to be affected by this problem. This claim is further 
supported by the calculation of the Schwingcr function 
on both coarse and fine lattices. We find that reflection 
positivity is clearly violated for both lattices; see, for 
example, Fig. 8. 

At large Euclidean time, C(t) is expected to reach zero 
from below. Unfortunately, this cannot be seen from our 
data shown in Fig. 7. Something similar has been ob- 
served independently using a different fermion formula- 
tion [28] . The reason might be the difference in the num- 
ber of lattice points in spatial and temporal directions. It 
is known [16, 17] that on an asymmetric lattice the gluon 
propagator can violate the infinite-volume continuum- 
limit tensor structure (Eq. (3)) in the (extended) time 
direction at very low momentum values. One effect of 
this is to alter the normalization of D(p = 0) relative to 
D(p > 0). As is customary, we assumed the continuum 
tensor structure and thus D(p = 0) has been normalized 
by N d — 4 and D(p > 0) by N d - 1 = 3. However, us- 
ing a slightly different normalization either for -D(O) or 
D(p > 0) the data in Fig. 7 could be easily shifted to- 
wards larger values such that (within errors) they equal 
zero at large t. Therefore, finite-volume effects could 
be responsible for C(t) not reaching zero from below at 
large t. 

Apart from C(t), it is also interesting to consider its ab- 
solute value and to compare our lattice data with corre- 
sponding DSE results by Alkofer et al. [35]. Those results 
have been reprinted on the right-hand side of Fig. 9. On 
the left-hand side we present our lattice result for \C(t)\ 
on the fine ensemble. The two figures are remarkably 



alike: not only does the zero crossing occur at almost 
the same place in Euclidean time, f~5 GeV -1 (about 
the size of a hadron), but the DSE results also correctly 
predict the small shift due to the inclusion of dynamical 
quarks. This shift is largest in our data for the smallest 
light quark mass. A similar result is obtained by using 
the coarse lattices. 



V. CONCLUSIONS 

We have extended a previous lattice study [24] of the 
Landau-gauge gluon propagator in full QCD to a finer 
lattice. The addition of quark loops has a clear, quanti- 
tative effect on the gluon propagator, but its basic fea- 
tures are unaltered. Good scaling behavior is observed 
for the gluon dressing function in both the quenched and 
unquenched cases with these improved actions. 

The violation of reflection positivity of the gluon prop- 
agator was investigated by calculating the real space 
propagator, or Schwinger function; the Landau gauge 
gluon propagator clearly violates positivity in both 
quenched and full QCD. 
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